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1. Introduction 



In the last decade, the market for non-volatile computer storage has experienced rapid growth. 
However, as the device size scales down, currently used NOR and NAND technology is 
facing limitations from tunnel oxide and electrostatic interactions between cells (Lacaita & 
Wouters, |2008|) . To solve the scaling problems, in the last fifteen years, a number of alternative 
Flash memory technologies have been proposed and studied. The phase-change memory, 
solid electrolyte memory, ferroelectric memory, magnetic memory and molecular memory 
using conducting molecules and carbon nanotubes all are promising can didates to replace 
technologies that are reaching their limit JChung et alll2010tlDraboidll2009l) . 

Vitreous materials involving chemical species from column VI other than oxygen are called 
Chalcogenide glasses. Chalcogenide glasses have long held interest for applications such 
as infrared detectors and optical fibers. Recently, chalcogenide glasses attracted further 
£^ ' attention due to their promising application in data storage devices. Two of the examples 
\ are phase-change memory materials based on tellurium alloys and solid electrolyte memory 
. materials based on metal-doped Ge-Se glasses. In this chapter, by using ab-initio molecular 
— i ' dynamics, we introduce the latest simulation results on two materials for flash memory 
devices: Ge2Sb2Tes and Ge-Se-Cu-Ag. Thi s chapter is a rev iew o f our previous work inclu ding 
. ■ some of our published figures and text in ICai et al.l J201(]|) and E rasai & Draboldl Em)) and 
also includes several new results. We organized the chapter as follows: we first introduce the 
, simulation m ethod used in a ll calculations in section|2j in section[3] we show the key findings 
on GeTSbT Te^Cai et al.l.l2010h: then we switch to the study on Ge-Se-Cu-Ag in sectionfll ( Prasai 
& DraboldT bOllI) . For both cases, after forming realistic atomistic models, we analyze the 
topology and electronic structure, predict their properties, and compare with experimental 
results. 

2. Simulation Method 



We use the ab-initio molecular dynamics method (ab-MD) to generate atomistic models. When 
such schemes are applied, the initial position of atoms are usually randomized and the system 
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is then annealed or relaxed to seek suitable local energy minima. The final model will be the 
one with a minimum total energy and agreeing with experimental measurements. One of 
the popular techniques based on ab-MD is the "Cook and Quench" method, for which the 
MD simulation is performed at a temperature well above melting point, which will force the 
system to lose memory of the initial configuration. Finally, the system is equilibrated at a 
lower temperature, like room temperature. Then an energy minimization is applied. Many 
realistic models are made by such a simple but powerful method. 

In our work, all of the calculations were carried out using periodic boundary conditions with 
the Vienna Ab-initio Sim ulation Package (VASP) . VASP is based on density functional theory 
using a plane wave basis jKresse & Furthmullerl.[l996|) . For Ge2Sb2Tes, we used the projector 
augment ed-wave (PAW) potentials and generalized gr adient approximation PBE (GGA-PBE) 
method fcresse & Toubertl, Il999l; iPerdew et ail, 1 19961) ; for Ge-Se-Cu-Ag, we used the local 
density approximation (LDA) for the exchange correlation energy in conjunction wi th the 
Vanderbilt Ultra Soft pseudopotentials dKresse & Hafnerl.ll994;|Perdew & Zungerl.ll98l|) . Both 
systems were annealed, equilibrated and cooled using molecular dynamic (MD) option of 
VASP and relaxation is carried out in conjugate gradient (CG) mode. Moreover, to obtain a 
better estimation for electronic gap, we applied Hartree-Fock (HF) calculation when analyzing 
the electronic structure of amorphous Ge2Sb2Tes. Though HF is known to exaggerate both 
the optical gap and charge fluctuation in the electron gas. These features are helpful to us for 
diagnosing the correlations between topology and electronic properties. 

3. Phase-Change Memory Material 

3.1 Background 

For Ge-Sb-Te (GST) alloys, there exists a rapid and reversible transition between crystalline 
and amorphous states. Controlled modification of electrical conductivity and optical 
properties of the transition is the ba sis for promising FLASH and optical memory devices. 
Akola and Jones dAkola & Tonesll200^) analyzed the structure of liquid and amorphous phases, 
and compared the electro nic structure with the crystal phase. In 2008, Hegedus and Elliott 
faegedus & Elliot! 120081) reproduced the crystal-amorphous transition by MD simulation, 
and they found that the rapid crystal growth was due to the presence of crystal fragments - 
four member square rings (so-called "seeds") in amorphous and liquid phases. Their work 
provided a way to track the dynamic chan ges of network top ology and electronic structure 
at the same time. Welnic and co-workers JWelnic et all |2007|) studied the origin of optical 
properties and argued that the optical contrast between amorphous and crystalline phases is 
due to a change in local order of Ge atoms. Despite this progress, the correlation between 
topology and electronic structure, most especially the origin of the change in the electronic 
gap, is still imperfectly understood. One of the challenges is the basic limitation of the LDA 
for estimating the gap. 

3.2 Model Preparation 

We began our work by creating amorphous Ge2Sb2Tes models by using the Vienna Ab-initio 
Simulation Packag e (VASP)- a plane-wave DFT code, using a PAW potential and the 
GGA-PBE method dHegedus & Elliot A 120081) . 63-atom amorphous Ge 2 Sb 2 Te 5 models with 
lattice constant 12.5 A were made as follows. The system was first melted and equilibrated at 
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Table 1. Mean coordinations, bond types and seeds(four member square rings) statistics at 500K. The 
result obtained at 1000K is listed in brackets (coordination cutoff=3.2A). 

Ntot N Te N Ge N sb 

Te 3.4(3.0) 20%(30%) 47%(41%) 33%(29%) 
Ge 4.6(4.3) 86%(71%) 5%(13%) 9%(16%) 
Sb 4.1(3.6) 69%(62%) 11%(20%) 20%(18%) 
N seed 18(1.8) 52%(10%) 69%(12%) 53%(10%) 

1000K, followed by a rapid quench to 500K with a quench rate of 16K/ps. Then the system 
was equilibrated for 20 ps and data collection began at lOps. For the crystal phase, 108-atom 
crystal Ge2Sb2Tes cells with lattice constant 21.316A are generated based on NaCl rock-salt 
structure with 10% vacancies: 60 Te atoms occupied the Na sites; 24 Ge atoms and 24 Sb 
atoms randomly occupied the CI sites which left 10 CI sites unoccupied. The system was then 
relaxed under zero-pressure till the minimum total energy was obtained. 

3.3 Result and Discussions 



3.3.1 Bond Statistics 




Distance (Ang) 

Fig. 1. Partial radial distribution functions for a-Ge2Sb2Te5. 

The calculated atomic coordinations for a-Ge2Sb2Te5(500K) a nd l-Ge^Sb^Te^lOOOK ) are listed 
in Table 1 with a 3.2A cut-off. These results are similar to lAkola & Tones! te007|) . although 
in our case, the mean coordination of Ge atoms is slightly increased after thermal quench 
and equilibration, which may be due to the higher equilibration temperature used and /or 
size artifacts for our smaller model. More highly-coordinated Ge and Sb (5-fold,6-fold) atoms 
appeared in the amorphous phase, which suggests that a near-octahedral structure may be 
formed (square-rings and 8-atom cubes). These results indicate that structural ordering is 
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enhanced in the amorphous phase relative to the liquid. Moreover, the number of wrong 
bonds (Te-Te,Ge-Ge,Sb-Sb and Ge-Sb) are decreased from 1000K to 500K which indicates that 
the chemical order is also improved. The average number of "seeds" (four member square 
rings) shows an increase in the amorphous phase and more than 50% of the atoms are involved 
in "seeds", compared to only 10% in the liquid phase. The calculated partial radial distribution 
functions are plotted in Fig. [10] The first peak in the Te-Ge and Te-Sb partials are located at 
2.81A and 2.92A. The shallow first minima imply that the coordination is sensitive to the 
cutoff value selected. The Te-Ge partial has a broad and weak second peak. However, the 
Te-Sb partial possesses a second peak with a maximum at 4.4A which indicates that Ge and 
Sb atoms differ in local environment relative to Te atoms. Regarding the homopolar bonds, 
there is a m ajor peak for the Sb-S b partial centered at 2.9A. The se results are simi lar to other 
simulations jAkola & Tones!. |2007|) and also experimental results JNatio et"ai],l201C]|) . 

3.3.2 Electronic Structure 



1 T 




Energy (eV) 

Fig. 2. Electronic densities of states projected onto different atomic species and orbitals. The Fermi level 
is at eV. 

The electronic structure is analyzed through the electronic density of states (EDOS) obtained 
from Hartree-Fock (HF) calculations. HF is used only to analyze the EDOS, not for forces 
and total energies. HF is known to exaggerate both the optical gap and charge fluctuations in 
the electron gas. These features are helpful to us here for diagnosing structural correlations. 
In the following discussion, the calculated EDOS is averaged over 1000 configurations from 
the last 2 ps when the cell is in thermal equilibrium at a fixed temperature of 500K. Finally, 
the averaged HF result of the amorp hous phase gives an e lectronic gap around 0.4eV which 
is wider tha n the DFT resul t-0.2eV dAkola & TonesL |2007|) and is closer to the experimental 
value-0.7eV jLee et al.l.l2005)l . Although the gap is still smaller than the experimental value, it 



Will-be-set-by-IN-TECH 



5 




Fig. 3. Projected EDOS on Ge atoms at tetrahedral and octahedral sites. "T" and "O" represent tetrahedral 
and octahedral sites. The "Ge-T/O site" plot only considers the contribution of Ge atoms to the EDOS, 
while the "Ge-T/O site with neighbors" plot contains the contribution of Ge atoms and their neighbors. 
The Fermi level is at eV. 

To correlate topology with electronic structure, we projected the EDOS onto different local 
sites and are able to attribute the electronic states to specific structural units. We first show 
the averaged EDOS for different species and orbitals in Fig. |2] The key findings are that, 
for all three species, p orbitals dominate the gap and tail states; if considering the species, 
Te-p,Sb-p,Ge-p,Ge-s and Sb-s are important in determining tail states and the magnitude of 
the gap (Fig. |2). To further correlate structural oddities with electronic states, we also sort 
atoms with specific features into different groups and accumulate the contribution to EDOS. 
We briefly report that groups forming homopolar or heteropolar bonds showed that there is 
a significant difference at a "deep gap" around -7eV below the Fermi level (OeV) in EDOS 
(atoms involved in heteropolar bonds form a bigger deep gap); however, atoms forming 
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homopolar bonds have a minor impact on tail states and the electronic gap near the Fermi 
level. Considering the coordination, for Te, 2-fold Te atoms contribute to a narrowed gap and 
conduction-band tail states appear; for Ge atoms, the contributions for 3, 4, 5 and 6-fold atoms 
are almost the same; for Sb atoms, the conduction-band tail of 6-fold Sb atoms is pushed to 
a low-energy level and the valence-band tail associated with 3-fold Sb atoms which satisfy 
the "8-N" rule is pushed into a higher energy region. While there are differences in electronic 
tail states and the gap value associated with coordination numbers, the influence is fairly 
weak. Similarly, sorting atoms involved in "seeds" or not also showed a minor impact on gap 
magnitude and tail states. 

We also considered the "umbrella flip" of Ge atoms. We compared the EDOS of Ge atoms 
sitting at octahedral sites (O-site) and tetrahedral sites (T-site), as we illustrate in Fig. [3] The 
projected EDOS on Ge atoms and their neighbors are all considered. The results indicate 
that 6-fold octahedral Ge and tetrahedral Ge have a similar local gap. However, 4-fold Ge at 
an octahedral site (4 neighbors with 90 degree angles) have both a shifted valence-band tail 
and conduction-band tail which may result in a bigger gap. Thus, from our result, sp 3 hybrids 
introduced by a Ge umbrella-flip may not be the reason for an increased gap in the amorphous 
phase, but the octahedral Ge existing in the amorphous phase at least would n ot increase the 
electronic gap. Analysis of GejSb2Te5 showed a similar result jRaty et aT1,l2010l) . 




Energy (eV) 



Fig. 4. AC Dielectric function of a-Ge2Sb2Tes. Due to finite size effect, the calculation can not predict 
valid value for small w (w < 2eV). 

Finally for this section, we show one last static property - the dielectric function e(to) in Fig.|4] 
The imaginary part of the dielectric function in three directions are plotted in Fig.|4] e reaches 
its peak for an energy of about 2.5eV and this spectrum is c omparable with both experimental 
and simulation result jRaty et alJ.l2uTollWuttig et al.l.l2007|) . Notice that due to the finite-size 



Will-be-set-by-IN-TECH 



7 



effect, the result is not valid for co — > 0. To obtain accurate results for small co, an extrapolation 
procedure is required. 



3.3.3 Dynamic Analysis 




Fig. 5. Dynamic change of temperature, gap value, number of wrong bonds and squares (seeds). 

Next, we performed a dynamic analysis for a-Ge2Sb2Tes. We tracked the structure and the 
electronic gap during a quench from 1000K to 500K with thermal equilibration at 500K (Fig.[3). 
Significant structural changes started to occur after 24ps (the temperature was then near 
640K). The number of homopolar bonds dropped, the number of 4-membered rings increased, 
an d the mean coordination increased. The changes in topology are similar to those reported 
by lHegedus & Elliottl J2008h and all these shifts signal an increase of both chemical order and 
structural order. The electronic gap, which we take to be the difference between LUMO and 
HOMO levels, increased overall, but we observed that there are considerable fluctuations, 
even for the well-equilibrated system. Local geometry may have huge consequences on the 

gap- 
To study how changes in the local environment at a specific site affected the electronic gap, we 
tracked a specific unit in the system during equilibration and we show such an evolution for 
both the topology and the electronic structure in Fig.[6]and Fig.[7] We mainly focused on one 
Ge atom which occupied a near-octahedral site (6 nearest Te neighbors with around 90 degree 
bond angles, indicated by a black arrow in Fig. [6) and its six nearest neighbors. We correlated 
their local bondings and electronic density of states for many time steps. Configurations 1 
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Config.4 



Fig. 6. Snapshots of topology changes for one Ge atom and its six neighbors (Ge-blue, Sb-brown, 
Te-green). The central Ge atom is identified by black arrows. The valence-band tail states appear in 
Config.2 and are localized on yellow atoms. The conduction-band tail states appear in Config.3 and are 
localized on black atoms. 




-0.08 0.67 
Energy (eV) 



Fig. 7. Instantaneous snapshots of EDOS correlated with the configurations of Fig. [5] A smaller gap 
appear in Config.2&3. The valence-band tail states (orange arrow) are associated with yellow atoms 
in Config.2 of Fig. [6] The conduction-band tail states (black arrow) are associated with black atoms in 
Config.3 of Fig.[6] The Fermi level is at eV. 
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and 4 exhibit the biggest gap. However, at intermediate steps between configurations 1 to 4, 
tail states appear. At configuration 2, a valence-band tail state was present and it was mainly 
localized on the central Ge atom and four of its nearest neighbors (yellow atoms in Config.2 of 
Fig. [6); at configuration 3, a conduction-band tail appears, mainly localized on the center Ge 
atom and two of its nearest neighbors (black atoms in Config.3 of Fig.[6). We should emphasize 
that from configurations 1 to 4, the whole network did not experience a major change, but the 
electronic gap fluctuates. Thus, the appearance of valence-band and conduction-band tails are 
strongly associated with distortions at this Ge site. Our simulations emphasize the dynamic 
nature of the electronic band tails in Ge2Sb2Te5. 

3.3.4 Relaxation Analysis for Crystal Phase of Ge^Sb^Je^, 




Fig. 8. Electronic density of states of crystal models projected onto different species of atoms. (Top-panel) 
Unrelaxed crystal model with vacancies. (Bottom-panel) Relaxed crystal model. 

In this section, we discuss relaxation effects for crystalline Ge2Sb2Tes with 10% vacancies. 
As mentioned above, the 108-atom cell was obtained based on the NaCl rock-salt structure. 
We show the electronic density of states for both unrelaxed and relaxed models in Fig. [8] 
obtained through HF calculations. For both models, Te atoms have a major effect on the 
valence tail which may be due to the vacancies; Sb atoms contribute more to conduction 
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Fig. 9. Change of LUMO, HOMO level, gap value and total energy during relaxation. 

tail. We could see clearly that the electronic gap opened up after relaxation. Moreover, we 
tracked the dynamic change of the Highest Occupied Molecular Orbital (HOMO) and Lowest 
Unoccupied Molecular Orbital (LUMO) and plot them in Fig. [9] It is clear from the plot that 
the total energy is reduced and both HOMO and LUMO levels are shifted. The HOMO level 
is pushed away by O.leV and the LUMO level is pushed up by 0.2eV. 

Since the crystal model has 10% vacancies, the relaxation actually introduced slight distortion 
into the network. The structural statistics indicate that the mean coordination of Te, Ge and 
Sb atoms all decreased. The mean coordination of Te are decreased from 4.8 to 4.28, Sb and Ge 
dropped from 6 to 5.47 and 5.23 correspondingly. The angle distribution, especially the X-Ge-X 
and X-Sb-X angle distributions, are also changed. This result indicates that the existence of 
vacancies and the distortion happened to the network will have a impact on gap. Thus, by 
controlling the concentration of vacancies and distortion, we may ob tained different elec tronic 
gap values. This result is similar to results on other Ge-Sb-Te alloys jWuttig et alj,l200^ . 

3.3.5 Conclusions on Ge^Sh^Ie^ 

We made Ge2Sb2Tes models with a 'quench from melt' method. HF calculations give a 
0.4eV electronic gap for the amorphous phase. We found that Te-p, Sb-p, Ge-p, Ge-s and 
Sb-s orbitals are most important to tail states. 6-fold octahedral Ge and 4-fold tetrahedral 
Ge give rise to similar gaps but 4-fold octahedral Ge results in a bigger gap with both 
shifted valence-band and conduction-band tails. The study also reveals a large fluctuation 
in gap value during thermal equilibration which is partially due to the appearance and 
disappearance of conduction-band and valence-band tail states. Such fluctuations could be 
associated with the local structural change/ distortion of Ge atoms, which introduce localized 
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tail states and have an impact on the electronic gap. Also, the relaxation analysis on crystal 
phase of Ge2Sb2Tes indicates that vacancies and distortions may play an important role in 
determining the electronic gap. 

4. Electrolyte Materials 
4.1 Background 

Electrolytes are materials with high ionic conductivity and high electrical resistivity. When 
doped with metals like Ag, chalcogenide glasses (e.g. Ge-Se) become solid electrolytes 
offering high ionic conductivities. Such electrolytes are getting attention for their 
technological import a nce w ith the application in "conducting bridge" (flash) memory devices 
jMitkova & KozickiL |2002|) . It has been believed that a variety of different coordination 
patterns of Ag + ions in the glassy host with tiny energy differences is the basic reason why 
Ag+ is mobile. Since the properties of chalcogenide glasses accrue from their structure, the 
knowledge of the structure of these glasses is an essential precursor for further study. From a 
material point of view it is interesting that an amorphous material should allow rapid motion 
of a transition metal ion through the network, and a great deal of energy has been devoted to 
understanding this phenomenon. Such diffusive processes in glasses have been studied for 
decades with a variety of experimental methods. There have also been several approaches to 
modeling such diffusive behavior. 

There have been a wide range of experimental studies on the atomic structure of the 
amorphous state of electrolyte material and some computer simulations, typically on Ge-Se 
glasses doped with transition metals. Ge-Se- Ag based electrolyt e materials have been stu died 
experim entally using various techniques. For example, X-ray {Piarristeguy et all l2000h and 
neutron jCuello et all |2007|; iDejus et all Il992h diffraction, and other experimental methods 
have been used to study the structure of Ge-Se-Ag glas s. There have al so been some 
computational studies to model the structure. Tafen et al. jTafen et all 12005^ reported two 
ab-initio models; (GeSe3)o.9Ag jand (GeSe3)o.85Ag .i5with short range order consistent with 
the experimental results. It has also been reported that Ag atoms prefer to sit at trapping 
center (TC) which is near the midpoint of a line joining two host atoms (Ge or Se) separated 
by a distance bet ween 4.7 and 5.2 A wit h the bond length of Ag to the host atoms ranging 
between 2.4-2.6 A jChaudhuri et alll2009|) for low Ag co ncentration. The simulat ion work has 
been also extended by introducing Cu into the network jPrasai & Drabold1,l201ll) . 

Beside structural studies, there have been quite a few studies on the conductivity of Ag doped 
chalcogenide glasses including both experimental and simulation work. Ag x (GeSe3)ioo-.v 
glasses have been particularly studied for the ioni c conductivity within a wide range of 
x (10 to 25%). Urena et al. iUrena et all l2005h predicted that the ionic conductivity 
follows an Arrhenius law. Tafen et al. presented a mole cular dynamics(MD) simulation on 
Ag. v (GeSe3)ioo-.Y with x = 10 and 15% iTafen et"aill2005|) at different temperatures. In recent 
work, we have also presented a MD simulations on these glasses with the addition of Cu and 
illustrated the motion of the ions on the accessible time scales (tens of picoseconds) (Prasai & 
Drabold. l201l|) . Some of the results will be discussed in the following sections. 
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Fig. 10. Comparison of total radial distribution functions and static structure factors for all amorphous 
models 

4.2 Simulation of properties of Electrolyte Materials 

The models of Ag- and Cu-doped chalcogenide glasses discussed here were generated using 
the melt-quenching method. A cubic supercell is constructed with a fixed volume and a 
fixed number of atoms in order to reproduce the experimental density according to the 
desired stoichiometry. The atoms were randomly placed in the supercell with minimum 
acceptable distance between two atoms set to 2A. The calculations were carried out under 
periodic bo undary condition using the Vienna Ab-initio Sim ulation Packaee(VASP)(Kresse & 
Furthmuller, Il996l) , with Vanderbilt ultrasoft pseudopotentials. We used the local density 



approximation (LDA) for the exchange correlat i on en ergy. The details of the model generation 
can be found in the reference lPrasai & Draboldl j201l|) . Beside the models discussed there, two 
more models (GeSe3)o.sCuo.2 and (GeSe3)o.8Cu .i Ag .i have been added to the discussion. 

4.3 Results and Discussion 



4.3.2 Structural properties 

Fig. [10] shows the calculated total radial distribution functions (RDFs) and structure 
factors for the models; g-(GeSe 3 ) .9Ag .i, g-(GeSe 3 )o.8Ag . 2 , g-(GeSe 3 ) .9Cu .iand 
g-(GeSe3)o.77Cuo.rj3Ago.2- The first peak of the RDF is the contribution from Ge-Se and 
Se-Se correlations whereas the second peak is due to Se-Se and Ge-Ag/ Cu correlationsfFig.[TT1 
and Fig.[12j. There is not much variation in the short range order (SRO) i.e. nearest neighbor 
distance and second nearest neighbor distance for the different models. We observed a 
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Fig. 11. Partial radial distribution functions for amorphous (GeSe3)o.9Ago.i (black) and 
(GeSe3)o.8Ag .2(red/ dashed line) 

slight change in the nearest neighbor distance for the Ag rich model and Cu rich model. The 
average bond length and the mean coordination numbers are presented in Table |2] We did 
not detect Ge- Ge bonds in any of our models as seen p reviousl y in g-(GeSe^1n oAgn i (T afen et 
al.,|200j). We also observed that both Ag and Cu preferred to have Se as neighbor with only 
16% of Cu/Ag bonded with Ge i n our models. These res ults are very close to bond lengths 
measured by Piarristeguy et al. jPiarristeguy et al.l.l2000|) . We also obtained the silver and 
copper coordination number for each model. The coordina tion number 3.1 of silver at 20% 
is as predicted(3.0) by Mitkova et al. jMitkova et al.l, [l999h . The coordination number 4.67 
of copper at 10% is m uch higher than 2.16 of silver (foun d to be 2.0 bv T afen et al. (Tafen et 
al., l2005l) ) for the same concentration. We detected a few 3-fold Ge and 3 and 4 fold Se that 
we interpret as a struct ural defect in our models. Detailed bond parameters can be found in 
IPrasai & Draboldli2uTll) . 

We also compared the static structure factors for our models (Fig.llOt. There is no significant 
change in the position of the first two peaks. We observed a weak peak in S(Q) slightly above 
1 A . This peak, which is a precursor to the first sharp diffraction peak (FSDP), varies 
as a function of Ag concentration an d the peak disappears as Ag concentration increases, 
also shown by Piarristeguy et al. JPiarristeguy et all 120031) . We did not observe any 
particular correlation contributing to this peak as the partial structure factors shows that 
the peak has contribution from all of the partials. We compared partial structure factors for 
(GeSe3)o.9 Ago i and (GeSe3)o.9Cug j and observed the only differences in correlation of Ag-Ag 
and Cu-Cu as well as in Se-Ag/ Cu. 
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Fig. 12. Partial radial distribution functions for amorphous (GeSe3)o.9Ago.i (black) and 
(GeSe3)o.9Cu .i (green/ thin line) 

Table 2. Short range order; nearest neighbor distance(NN), next nearest neighbor distance(NNN) and 
mean coordination number(CN). 



NN(A) NNN(A) CN 



(GeSe 3 ) .9Ag al 


2.49 


3.75 


2.50 


(GeSe 3 )o.8Ag .2 


2.51 


3.80 


2.92 


(GeSe 3 ) .77Cu .o3Ago.2 


2.45 


3.80 


2.9 


(GeSe 3 ) .9Cu al 


2.40 


3.83 


2.8 



We performed thermal MD simulation at 1000K for 25ps in order to obtain well-equilibrated 
liquid systems. We calculated the total and partial radial distribution functions (RDF). 
The RDFs are averaged over the last 2.5 ps. The major peak positions in total RDF are 
2.45 A for (GeSe 3 ) .9Cu . 1/ 2.48 A for (GeSe 3 ) . 9 Ag 1 and 2.53 A for (GeSe 3 )o. 8 Ag a2 and 
(GeSe 3 )o.77Cuo.o3Ago.2- We present partial radial distribution functions in Fig. [13] 
showing Ge-Ge, Ge-Se, Se-Se and Se-Ag/Cu correlations. All of our models except 
(GeSe 3 )o.9Ago.i(2.6A) confirm the presence of Ge-Ge homopolar bonds with peak position 
at 2.71 A in contrast with the glass. We also observed Se-Se and Ge-Se bond distances of 2.47 A 
and 2.50 A, respectively. We observe no concentration dependence on the first peak position 
of Ge-Se,Se-Se and Se-Ag/ Cu correlations. The major contribution to the first peak of the total 
RDF is from Ge-Se,Se-Se and Se-Ag/Cu correlations with Se-Ag/Cu correlation causing the 
shifts on the first peak positions. The second peak of the total RDF is mainly due to Se-Se 
correlation. 
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Fig. 13. Comparison of partial radial distribution functions for all liquid models at 1000K 
4.3.2 Ion dynamics 

We studied the dynamics of Ag and Cu ions in the GeSe3 host by computing the mean square 
displacement (MSD) for each atomic constituent as: 

i N„ 

(r 2 (t)) a = ^E(\m-m\ 2 ) a) 

where the quantity in {) is the calculated statistical average over the particular atomic species 
a. We carried out constant temperature MD calculations at three different temperatures 300K, 
700K and 1000K in order to study ion dynamics in our the amorphous as well as the liquid 
systems. 

4.3.2.1 Amorphous Ge-Se-Cu-Ag 

As expected, at 300K none of the ions showed measurable diffusion. In order to investigate 
the diffusion in the solid place, we chose T = 700K and present the MSD for each species 
for each system calculated at this temperatu re in Fig. [HI At 70 0K Ag+ ions show significant 
diffusion consistent with the previous result flafen et al.l . 12005ft in contrast to Cu ions that do 
not diffuse much. To elucidate the diffusion of these ions we examine the trajectories for 
20ps. Fig. [15] and [16] show two dimensional projections of the trajectories of the most and 
the least diffusive ions in (GeSe3)o.9Ago.i and (GeSe3)o.9Cuo.i. The trajectories illustrate the 
wide range of diffusion for the ions with displacement ranging 1A-3.87A in (GeSe3)o.9Ag .i, 
2A-6.71A in (GeSe3)o.sAgo.2 and 1A-3.74A in (GeSe3)o.9Cu j. For the mixed-ion model 
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Time (ps) 

Fig. 14. Mean square displacement of atoms in amorphous (GeSe3)o.9Ag .i, (GeSe3)o.sAgo.2/ 
(GeSe3)o.77Cuo.o3Ag .2and (GeSe3)o.9Cu .i(top to bottom respectively) glasses at T = 700K. Ag(black) 
Ge(green), Se(red) and Cu(blue) 

(GeSe3)o.77Cuo.o3Ago.2, this displacement ranges between 1.73A-2.82A for Cu and 1.41A - 
8.06A for Ag. For Ag rich models more than 60% of the ions exhibit displacements greater than 
the average displacement (2.36A in (GeSe3)o.9Ag .i and 4.47A in (GeSe3)o.sAgo.2) whereas 
for Cu, the majority has displacement smaller than the average(2.11A). The wide range of 
diffusion can be attributed to variation in the local environment of the ions. To illustrate this 
we calculated the local densities of the most and the least mobile ions. We employed a sphere 
of radius 5. OA around the ion and calculated the mean density of atoms inside the sphere. We 
observed that the most diffusive ion is located in the region with lower local density. In other 
words the most mobile ions have the wider variation of the local density as compared to that 
of the least mobile ion. 

4.3.2.2 Liquid Ge-Se-Cu-Ag 

One of the essential properties of a liquid is the high diffusivity of atoms in the system. To 
illustrate this, we calculated the mean square displacements for each species at 1000K in all of 
our models. The diffusion plots as presented in Fig. [17] shows that the MSD of each species 
increases rapidly as compared to that at 700K. We observe Ag diffusion still significantly larger 
than the host particles however; Ge and Se atoms are also diffusing rapidly. As before Cu still 
does not show high diffusion as Ag does compared to the host atoms. 
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Displacement (A) 



Fig. 15. Trajectories of the most and the least diffusive Ag ions at 700K as a function of time in amorphous 
(GeSe 3 )o. 9 Ago.i. 



Based on the plots we calculated diffusion coefficients using Einstein relation <Chandlerl.[l987h . 
The Einstein relation for self-diffusion is given by: 



(\r i (t)-r l (0)\ 2 }=6Dt + C 



(2) 



where C is a constant and D is the self-diffusion coefficient. The conductivity can be calculated 
from the equation 

ne 2 D 

where n is the number density of ions. The temperature dependence of the diffusion is shown 
in Fig. [18] and the values of diffusion coefficients and conductivities at different temperatures 
are presented in Table [3] We did not find experimental results for the conductivity o f Cu ions; 
however Ag conductivity is close to ones reported by Urena et al. jUrena et all 120051) . 

4.3.3 Trap centers and hopping of ions 

To illustrate the different ionic transport properties of Ag and Cu, it is essential to study the 
local environment of Ag and Cu in our models. Fig. [19] shows the local environment for Ag 
and Cu in (GeSe3)o.9Ag .i and (GeSe3)o.9Cu j respectively. In the relaxed networks, most of 
the Ag ions(58.3%) are found to occ upy the trap centers, between two of the h ost sites as also 
predicted by the previous workers jChaudhuri et all 120091: fTafen et all Eoolfr but this is not 
the same case with Cu. Cu is always surrounded by more than two host atoms that makes 
the traps for Cu more rigid than for Ag. In Ag rich systems at 300K, we observed that Ag is 
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most mobile Cu 

least mobile Cu 




a 

Displacement (A) 

Fig. 16. Trajectories of the most and the least diffusive Cu ions at 700K as a function of time in amorphous 
(GeSe 3 )o.9Cuo.i. 

Table 3. Self diffusion coefficient D and conductivity a at 300K, 700K and 

1000K for (GeSe 3 )o.9Ag . 1 (10%Ag) / (GeSe 3 ) a8 Ag a2 (20%Ag), (GeSe 3 ) .9Cu al (10%Cu) and 
(GeSe 3 )o.77Cuo.o 3 Ago.2(0.77%Cu) 



T(K) D(cm 2 /s) ^(Scirr 1 ) 

This work Expt.Urena et al. (2005) 



10%Ag 


300K 


1.15xl0- y 




2.63x10 


-5 


1.3xl0- b 




700K 


4.53 xlO" 6 




4.44x10 


-2 


2.07X10- 2 




1000K 


1.23 xlO" 5 




8.45x10 


-2 


8.98xl0- 2 


20% 


300K 


1.16xl0- 8 




5.3x10- 


4 


7.5xl0- 5 




700K 


1.20 xlO" 5 




2.35x10 


-1 


6.57x10-2 




1000K 


2.53xl0- 5 




3.47x10 


-1 


2.584 xlO- 1 


10%Cu 


300K 


7.3xl0- 10 




1.67x10 


-5 






700K 


3.3 xHT 6 




3.23x10 


-2 






1000K 


1.13xl0" 5 




7.75x10 


-2 




0.77%Cu 


300K 


D Ag =1.06xlO~ 


-8 


4.85x10 


-4 








D c „=7.16xKT 


- c > 


1.63x10 


-5 






700K 


D Air =1.30xlO^ 


-5 


2.54x10 


-1 








D a ,=1.16xKT 


-6 


3.8x10- 


3 






1000K D A „=2.42xlO~ 


-5 


3.32x10 


-1 








D Cu =5.24xl0- 


-6 


1.2x10- 


2 
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Fig. 17. Mean square displacement of atoms in liquid (GeSe3)o.9Ago.i, (GeSe3)o.sAgo.2/ 
(GeSe3)o.77Cu .o3Ago.2and (GeSe3)o.9Cu .i(top to bottom respectively) glasses at T = 1000K. Ag(black) 
Ge(green), Se(red) and Cu(blue) 

basically trapped with only a few hopping events. At 700K the lifetime of the trap decreases 
and hopping occurs. We observed the lifetime of the traps varying from lps - 3.5ps. However 
at 1000K we failed to observe well defined hopping events because of the high the diffusion of 
the host itself. In the Cu rich system the story is completely different. Even at 700K we could 
observe only a few hopping events with much larger trap life time. It has also been shown 
by previous workers that the nature of trap or cage depends mainly on coordination number, 
nearest neighboring distance and angular distribution of the nearest neighbors (Kraemer & 
Naumis, l2008l) . The low coordination number of Ag makes it easy to escape the trap whereas 
for Cu, high coordination number, smaller neighbor distance and a more uniform angular 
distribution makes it more difficult to escape from the trap. 



4.3.4 Mixed Ion Conductivity 



One big challenge in these materials is to fully understand the effect on the dynamic properties 
such as ionic conductivity when one of the mobile ion is partially substituted by another type 
of mobile ion. There is a non-linear change in ionic mobility when two or more than two 
types of mobile ions are mixed in ion conducting glasses and crystals, and the effect is known 
as mixed ion effect. This section reveals that the mixed ion effect in Ag and Cu doped GeSe3 
glasses is present in our simulation. Constant temperature MD simulations were carried out 
in (GeSe3)o.8(Ag 1 _j.Cu v )o.2 where x = 0,0.5 and 1 at two different temperatures of 700K and 
1000K. The calculated ion conductivities are presented in Fig. [20] The figure shows a drastic 
drop in the ionic conductivity when both Ag and Cu ions are present in the system. This 
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Time (ps) 

Fig. 18. Temperature dependence of conductivity of ions for different models) 




Fig. 19. Local environments of Ag atoms(top) and Cu atoms (bottom). Black, green,blue and yellow 
colored atoms respectively represent Ag, Se, Ge and Cu) 
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Fig. 20. Mixed ion effect: comparison of ion conductivities (GeSe3)o.8(Agi_ x Cu I )o,2 glasses as a function 
of x. 

result implies a mixed ion effect in Ag/Cu doped chalcogenide glass, where Ag + conduction 
is greatly reduced by the presence of Cu + . It is encouraging to see a mixed-ion effect in our 
simulations; its atomistic origin is under study. 

4.4 Conclusion: Fast ion conducting glasses 

We prepared different Ag and Cu doped GeSe3 glass and liquid models by ab initio simulation 
using the 'melt-quench' method and analyzed their structural and electronic properties. We 
also simulated dynamics of Ag and Cu ions using molecular dynamics. We were able to 
reproduce structural data as provided by X-ray diffraction. From the electronic density of state 
we observed that the increase in Ag concentration widens the optical gap whereas increase in 
Cu concentration narrows the gap. We were also able to see the metallic behavior for the 
liquid systems with the gap closing completely at 1000K. We were able to show the diffusion 
of the ions even in our time scale and predict the conductivity close to the experimental data. 
We also studied the trap and found that Cu traps are more rigid that those for Ag making very 
hard for Cu to diffuse. 

5. Conclusion 

By using molecular dynamic simulations, we generate the atomistic models of Ge2Sb2Tes and 
Ge-Se-Ag-Cu and analyzed their topology and electronic structures. Both phase-change 
and electrolyte solid materials show the promising properties as candidates to replace the 
contemporary technologies in Flash memory. With further development, we believe the new 
generation of the computer storage device will eventually appear with much smaller size, 
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higher speed and more reliable features. We show that computer simulation can lend insight 
into promising technologies. 
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